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Abstract: VERITAS is a system of four imaging Cherenkov telescopes located at the Fred Lawrence 
Whipple Observatory in southern Arizona. We present here results of detailed Monte Carlo simulations 
of the array response to extensive air showers. Cherenkov image and shower parameter distributions are 
calculated and show good agreement with distributions obtained from observations of background cosmic 
rays and high-energy 7-rays. Cosmic-ray and 7-ray rates are accurately predicted by the simulations. The 
energy threshold of the 3-telescope system is about 150 GeV after 7-hadron separation cuts; the detection 
rate after 7-selection cuts for the Crab Nebula is 7.5 7s/min. The three-telescope system is able to detect 
a source with a flux equivalent to 10% of the Crab Nebula flux in 1.2 h of observations (5 a detection). 



Introduction 

The complete VERITAS array of four imaging 
Cherenkov telescope has been in operation since 
early 2007 at the Fred Lawrence Whipple Obser- 
vatory in southern Arizona. Several TeV 7-ray 
sources have already been detected and extensive 
reports about these observations can be found at 
this conference (see [1] and references therein). 
A detailed description of the first VERITAS tele- 
scope can be found in [2], and a full description 
of the VERITAS project in [3]. Results regarding 
a comparison of single telescope parameters from 
simulations and data can be found in [4]. This pa- 
per concentrates on the performance of the array of 
telescopes. 

Monte Carlo simulations 

The VERITAS collaboration uses several different 
simulation chains in order to be independent of 
possible systematic effects arising from individual 
simulation packages. The shower generators gen- 
erally used are KASCADE [5] and CORSIKA [6], 
while the response of the array of telescopes is sim- 
ulated with three different packages (ChiLa, KAS- 
CADE, GrISUDet [7]). A description of the results 
from all of the simulation packages is beyond the 
scope of this paper, and so we concentrate on the 
CORSIKA/GrlSUDet chain. 



CORSIKA version 6.501 is used with the hadronic 
interaction models QGSJet, for primary energies 
above 500 GeV, and FLUKA at lower energies. 
Gamma-ray, hydrogen and helium nuclei induced 
air showers are simulated in an energy range from 
50 GeV (30 GeV for hydrogen) to 30 TeV at differ- 
ent elevations. Spectral indices are taken for cos- 
mic rays from fits to balloon measurements [9, 10]. 
For 7-rays a distribution similar to the energy spec- 
trum of the Crab Nebula has been simulated. The 
shower cores are scattered randomly on a circu- 
lar area with a radius of 600 m around the center 
of the telescope array. The isotropic distribution 
of the cosmic -ray incident angles is simulated by 
randomizing the shower directions in a cone of ra- 
dius 3.5° around the pointing direction of the tele- 
scopes. 

Measurements of the atmospheric properties at the 
site of VERITAS are in progress [8]. The calcu- 
lations described here use the U.S. standard atmo- 
sphere, which does not always reflect the proper- 
ties of the atmosphere in southern Arizona. The 
photoelectron rate per PMT is measured to be 
100-200 MHz, corresponding to a night sky back- 
ground rate of 2.8 x 10 12 photons • m~ 2 s _1 sr _1 , 
which is used in the simulations. 

The telescope and array simulations consist of 
three parts: the propagation of Cherenkov pho- 
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tons through the optical system, the response of 
the camera and electronics, and the local and array 
trigger system. The geometrical properties of the 
optical system are fully implemented; misalign- 
ment of the mirrors and their surface roughness 
are taken into account. The camera for each tele- 
scope consists of 499 pixels with 1 1/8" phototubes 
and light cones. The response of the PMTs to sin- 
gle photons has been measured; the single photo- 
electron pulse has a rise time of 3.3 ns and a width 
of 6.5 ns [2]. In the simulations a signal in a PMT 
is created by summing single photo-electron pulses 
with appropriate amplitude and time jitter applied. 
Electronic noise and all efficiencies, including mir- 
ror reflectivities, geometrical, quantum, and col- 
lection efficiencies, and losses due to signal trans- 
mission have been modeled. The pulses are digi- 
tized with 2 ns sampling and with a trace length of 
24 samples, reflecting the properties of the FADC 
system. The trigger simulation utilizes a simpli- 
fied model of the constant fraction discriminator, 
which is the first stage of the VERITAS multi-level 
trigger, and a full realization of the pattern trig- 
ger, requiring three adjacent pixels above thresh- 
old in a time window of 5 ns. The currently used 
trigger threshold of 50 mV corresponds to about 
4.8 photoelectrons. The third level of the trig- 
ger system, the array trigger, activates the read 
out of the telescopes when at least two telescopes 
pass the pattern trigger requirements in a time win- 
dow of 100 ns. The analysis steps, which include 
pedestal calculation, image cleaning, image pa- 
rameterisation, source reconstruction, and calcu- 
lation of mean scale variables for 7/hadron sepa- 
ration are exactly the same for simulated and real 
data and are described elsewhere [11, 12]. 

Data 

For the following comparison of data with Monte 
Carlo simulations, 7-ray candidates are extracted 
from five hours of observations of the Crab Neb- 
ula in January and February 2007. The data were 
taken with an array of three telescopes (Telescopes 
1,2, and 3) in wobble mode (0.5° to 1.4° offsets) at 
a telescope elevation of about 70° in good weather 
conditions. 7-ray candidates were selected from 
the data using cuts on the shape of the shower im- 
ages (so-called mean scaled width and length cuts, 
see [11]) and on the shower direction (0 2 , the 
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Figure 1: Comparison of integrated size in digital 
counts of the third highest pixel in images in Tele- 
scope 3 for simulations (open symbols) and obser- 
vation data (closed symbols). 
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Table 1: Comparison of rates of 7-ray candidates 
after cuts for different telescope configurations and 
wobble offsets between data and simulations for 
observations of the Crab Nebula at elevations of 
about 70°. The errors for the rate derived from 
real data is statistical only; the stated errors for 
the Monte Carlo rate reflect the uncertainty of the 
spectral parameters of the assumed Crab Nebula 
energy spectrum. 

squared difference between reconstructed shower 
direction and source position in the sky). These 
cuts result in a 52 a detection of the Crab Nebula 
and provide ~ 1500 7-ray candidate events. 

Comparison of Data with Monte Carlo 
simulation 

The pattern trigger of VERITAS requires at least 
three adjacent pixels above threshold in a certain 
time window. The pixel with the third highest sig- 
nal size in a pattern therefore determines the acti- 
vation of the readout of the telescope. The com- 
parison between data and simulations of the dis- 
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Figure 2: Comparison of mean scaled width for 
7-ray and background data from simulations and 
observations. The distributions have been scaled 
to similar fluxes. The grey area indicates the range 
of values allowed by the 7-hadron separation cuts 
(same in Figure 3). 



Figure 3: Comparison of 2 distribution from sim- 
ulations (open symbols) and real data (closed sym- 
bols). The lines show fits to each of the distribu- 
tions (dashed line for simulations, continuous line 
for data). The distributions have been scaled to 
similar fluxes. 



tribution of the integrated charge measured in this 
pixel over many showers (Figure 1) demonstrates 
that the second level trigger simulation reflects the 
real system accurately. The event rate after basic 
quality cuts (requiring the successful reconstruc- 
tion of shower direction and impact parameter) is 
22.3±0.3 Hz for data and 20.6±0.4 Hz for simu- 
lated cosmic -ray showers. The typical dead time 
of the 3-telescope system of 6-7% is taken into ac- 
count. Before comparing these numbers for 7-rays, 
the parameters used for the suppression of back- 
ground events have to be tested. Figure 2 shows the 
mean scaled width distributions for measured and 
simulated 7-ray and background events, Figure 3 
the 2 distributions for data and simulations. The 
good agreement in both figures indicates that the 
simulations reflect the real system well. The rate of 
detected 7-rays from the Crab Nebula depends on 
the applied cuts. Those applied here are relatively 
hard: about 50% of all 7-rays are accepted and 
more than 99.9% of all cosmic -rays are rejected. 
Table 1 lists 7-ray rates for different array config- 
urations and wobble offsets for measurements and 
simulations. The off-axis acceptance of the sys- 
tem is particularly important for the detection of 
extended sources or during sky surveys. Figure 4 
shows that the acceptance is above 80%, relative 
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Figure 4: Relative 7-ray rate as function of dis- 
tance to camera center for data (black symbols) and 
simulations (red symbols) and 7-ray like showers 
(blue crosses). 



to a source at the center of the camera, for offsets 
from the camera center of up to 0.7° and is well 
described by the simulations. 

The energies of the primary 7-rays are recon- 
structed using lookup tables [12]. The bias in the 
reconstructed energy is, as Figure 5 shows, well 
below 10% for energies above 160 GeV. The effec- 
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Figure 5: Bias in energy reconstruction vs true en- 
ergy Emc for three different elevations. The error 
bars show the width of corresponding bias distri- 
bution (for the array configuration T1+T2+T3 and 
an observation mode of wobble offset 0.5°). The 
dotted lines indicate a bias of 10%. 
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Fi gure 6: Effective area vs true energy E^frj for 
different elevations (for the array configuration 
T1+T2+T3 and an observation mode of wobble 
offset 0.5°). 

tive area of the system is uniform above these en- 
ergies (Figure 6) and extends well beyond 10 TeV. 
The current analysis energy threshold of the sys- 
tem, defined as the position of the peak of the en- 
ergy spectrum of the source convolved with the ef- 
fective area curve of the detector, is about 150 GeV 
for observations near the zenith. It increases to 300 



GeV at 40° zenith angle. The three-telescope sys- 
tem is able to detect a source with a flux equivalent 
to 10% of the flux of the Crab Nebula in 1.2 h of 
observations (5 a detection). 

Conclusions 

Monte Carlo simulations of the VERITAS array 
of Cherenkov telescopes show that the system is 
well understood and accurately modeled. The 
agreement between the predicted and actual per- 
formance of the system, as well as the detection of 
several sources of high-energy 7-rays during the 
construction phase of the project, demonstrates the 
high potential of VERITAS for making detailed 
morphological and spectral studies of astrophysi- 
cal 7-ray sources from the northern hemisphere. 
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